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I.  INTRODUCTION 


To  properly  describe  chemical  reactions 
in  non-permeable  porous  media,  the  "pore  tree" 
was  introduced  by  Simons  and  Finson  (1979) 
and  Simons  (1982).  This  pore  structure  was 
developed  via  analogy  to  the  kinetic  theory  of 
gases:  the  pore  length  is  analogous  to  the  mean 
free  path  and  under  the  assumption  that  the 
pore  aspect  ratio  (pore  length  to  radius)  is  a 
constant,  a  pore  size  distribution  was  obtained 
that  has  been  confirmed  for  coal,  coal  char, 
sorbents,  catalysts  and  even  for  kidney  stones 
from  both  men  (low  porosity  oxalate)  and 
women  (high  porosity  phosphate).  The  pore 
tree  was  statistically  derived  from  the  pore  size 
distribution  and  allows  the  orderly  migration  of 
a  reactant  gas  from  the  large  pores  to  the  small 
pores  (see  Fig.  1).  A  detailed  description  of 
the  pore  tree  and  the  coupled  transport  and 
chemistry  is  given  by  Simons  (1982,  1983a). 

The  spatially  dependent  transport/reaction  equations  are  solved  for  a  single  pore  tree  and  then 
the  total  contribution  of  all  trees  (of  all  sizes)  in  the  system  is  obtained  by  summing  the 
contribution  of  each  tree  that  reaches  the  exterior  of  the  system.  This  is  distinct  from  the 
"bulk"  transport  equations  that  integrate  over  all  pores  at  a  fixed  point  in  space  before  they 
integrate  spatially.  The  "bulk"  transport  approach  inherently  assumes  instantaneous  mixing 
between  large  and  small  pores  at  a  single  point  in  space  (Random  Pore  Models)  while  the 
pore  tree  relaxes  this  assumption  with  a  structure  that  is  developed  statistically  and  is 
relatively  easy  to  work  with  mathematically. 

In  all  previous  applications  (sorbents,  catalysts  and  coal  combustion)  the  pore  tree 
represents  an  isolated  sub-structure,  allowing  diffusion  into  and  out  of  the  porous  media 
without  permitting  convection  through  the  media.  Under  the  current  ARO  program,  Simons 
Research  Associates  is  adapting  its  pore  tree  model  to  describe  the  permeable  pore  structure 
and  convective  transport  characteristic  of  the  subsurface  flow  of  water  in  soil,  the  dispersion 
of  contaminants  and  the  in-situ  remediation  of  contaminated  sites.  The  random  nature  of  the 
pore  structure,  which  formed  the  basis  of  the  statistical  derivation  of  the  pore  tree,  is  applied 
to  the  underground  structure  of  porous  soil,  sand  and  rocks.  The  statistical  description  of  the 
pore  tree  is  extended  to  describe  the  concept  of  permeability,  thus  allowing  convective  flow 
as  a  transport  mechanism  in  addition  to  the  diffusive  transport  already  considered.  Physically, 
this  amounts  to  statistically  determining  the  "branches"  that  are  common  to  several  trees  to 
allow  percolation  through  the  large  scale  (mobile)  structure  in  addition  to  diffusion  through 
the  smaller  scale  (immobile)  structure.  While  permeability  is  dominated  by  the  largest  pores, 
it  is  important  to  determine  the  level  of  convection  that  is  occurring  at  the  intermediate  scales 
in  order  to  accurately  describe  transport.  The  resulting  pore  structure  will  provide  an  analytic 
description  of  an  underground  network  upon  which  transport  and  coupled  chemical  reactions 
may  be  accurately  superimposed. 


rp  =  Appropriate 
Mean  Radius  at  x 


Particle  Surface 


Figure  1.  The  Pore  Tree 
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II.  ISOLATED  PORE  TREE:  THE  STRUCTURE 


Following  the  pore  structure  theory  of  Simons  and  Finson  (1979)  and  Simons  (1982), 
consider  a  spherical  porous  particle  of  radius  a,  containing  pores  of  length  lp  and  radius  rp. 
The  pore  dimensions  range  from  a  microscale  of  the  order  of  Angstroms  to  a  macroscale 
which  is  a  significant  fraction  of  the  particle  radius.  The  radius  of  the  largest  pore  is  denoted 
by  rmax  and  is  given  by 


^=2fl61/3/3  K0  (1) 

where  0  is  the  total  porosity  of  the  particle  and  K0  is  a  constant  of  integration,  approximately 
equal  to  five,  which  relates  the  pore  length  to  its  radius 

V~K>rJ*'p  <2> 


The  radius  of  the  smallest  pore  is  denoted  by  rmin  and  is  given  by 


9fr  (3) 

where  ps  is  the  density  of  the  solid  matrix,  sp  is  the  specific  internal  surface  area  (several 
hundred  m2/g),  and 

P=ta(r„AJ  (4) 


The  particle  contains  a  continuous  distribution  of  pore  sizes  from  rmin  to  rmax  .  The 
number  of  pores  within  an  arbitrary  plane  of  cross-sectional  area  A  and  with  radius  between 
rp  and  rp  +drp  is  denoted  by  g{r^Adrp.  The  pore  distribution  function  g{r^  is  given  by 

J(rp=0/2^prj 

where  g(r p  indicates  an  average  over  all  inclination  angles  between  the  axis  of  the  pore  and 
the  normal  to  the  plane.  Due  to  the  random  orientation  of  the  pores,  the  intersection  of  a 
circular  cylinder  with  a  plane  is  an  ellipse  of  average  area  27trp2.  Hence,  the  porosity  is  the 
27trp2  moment  of  g{r^  and  the  internal  surface  area  is  the  47trp  moment  of  g(r^ .  The 
expression  for  g(rp)  was  derived  (Simons  and  Finson,  1979)  from  statistical  arguments  and  has 
been  validated  through  extensive  comparison  of  the  predicted  volume  and  surface  area 
distributions  with  mercury  intrusion  data  (Stacy  and  Walker,  1972).  This  has  been 
accomplished  for  coal,  char  derived  from  that  coal  (Kothandaraman  et.al.,  1984),  sorbents, 
catalysts  and  even  kidney  stones  from  both  men  (oxalate)  and  women  (phosphate). 


3 


A  characteristic  feature  of  the  l/rp3  distribution  depicts  that  the  pore  volume  between 
r^  and  rp  increases  linearly  with  the  natural  log  of  rp.  It  is  the  functional  form  of  this 
relationship, 


Pore  Volume  « 


/  rpS(rj)drp  «ln rp 


rmin 


(6) 


that  depicts  the  inverse  cubic  dependence  of  g(r p  on  rp.  A  linear  display  of  mercury  intrusion 
volume  vs.  ln(rp)  always  infers  a  l/rp3  distribution. 

The  number  of  pores  within  the  bulk  volume  V  whose  pore  radius  is  between  rp  and 
rp+drp  may  be  defined  by  Vj{rp)drp .  The  pore  volume  is  expressed  as  the  Jtrp2lp  moment  of 
/(rp  and  the  internal  surface  area  is  the  27trplp  moment  of  /(rp .  The  pore  size  distribution 
functions  (  /(rp  and  g(r p  )  are  clearly  not  independent.  The  definitions  of  porosity  and 
internal  surface  area  infer  that  /(rp  is  related  to  g(r p  by 

(7) 

Equation  (7)  simply  states  that  the  probable  number  of  pores  intersecting  an  arbitrary  plane 
increases  with  the  length  of  the  pore  and  with  the  density  of  pores. 

The  length  of  a  pore  is  determined  by  an  arbitrary  intersection  with  another  pore  and 
is  expressed  (Simons  and  Finson,  1979)  as  a  collision  integral  over  the  pore  distribution 
functions.  The  analysis  suggests  that  ,  g(r p  and  /(rp  are  proportional  to  rp  ,  l/rp3  and  l/rp4 
respectively.  The  constants  of  proportionality  are  obtained  from  integral  constraints,  i.e.,  the 
total  porosity  and  internal  surface  area  contained  in  the  pore  structure.  The  expression  for 

is  given  by 


*P  K/p 


where  the  constants  were  defined  above. 

The  pore  volume  distribution  corresponding  to  these  distribution  functions  is  similar  to 
that  utilized  in  the  random  pore  model  (Gavalas,  1980  &  1981).  However,  the  pore  tree 
model  and  the  random  pore  model  differ  dramatically  in  their  choice  of  the  pore  aspect  ratio 
(length  to  diameter)  and  its  implications  with  respect  to  pore  branching.  The  random  pore 
model  allows  a  single  pore  to  connect  two  larger  pores.  This  picture  lends  itself  to  the 
idealization  of  instantaneous  mixing  between  the  pores  and  requires  that  the  pore  aspect  ratio 
be  of  the  order  of  one  hundred.  The  pore  tree  theory  uses  data  for  r^  to  imply  (via  K0)  that 
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all  pores  possess  an  aspect  ratio  of  the  order  of  ten.  Hence,  small  pores  may  connect  to  larger 
pores  only  on  one  end  and  all  pores  must  branch  from  successively  larger  pores  like  a  tree  or 
river  system. 

Each  pore  that  reaches  the  exterior  surface  of  the  particle  is  depicted  as  the  trunk  of  a 
tree.  The  size  distribution  of  tree  trunks  on  the  exterior  surface  of  the  particle  is  denoted  by 

g(rt)4na2drt  where  g(r)  is  functionally  identical  to  g(r^.  Each  trunk  of  radius  rt  is 
associated  with  a  specific  tree-like  structure.  Let  Nt  be  defined  as  the  branch  distribution 
function  where  Ntdrp  is  the  number  of  pores  of  radius  rp  (within  size  range  drp)  in  a  tree 
whose  trunk  radius  is  rt .  The  total  number  of  pores  of  radius  rp  in  a  sphere  of  radius  a 

may  be  expressed  as  4/3 n  #  r p)drp  or>  as  the  sum  °f  all  pores  of  radius  rp  contained  within 
every  tree  in  the  porous  sample,  plus  all  pores  of  radius  rp  that  are  themselves  the  trunk  of  a 
tree.  Hence, 


4 

—na 
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J  Ng(r)4%a2drt  +  4rca2g(rp) 


(9) 


where  g(r)  is  the  number  of  tree  trunks  per  unit  external  area  of  the  porous  sample  and  only 
those  trees  whose  trunk  radius  is  greater  than  rp  may  contain  a  pore  of  radius  rp .  Using  the 
previously  derived  expressions  for  r^,  g{r^  and  f{r p,  Eq.(9)  is  identically  satisfied  by 


»r  3,  4 

Nt=rjrp 


(10) 


The  branch  distribution  function  completely  characterizes  the  pore  tree.  The  internal 
surface  area  and  pore  volume  associated  with  each  pore  tree  are  denoted  by  St(r,)  and  Vt(rt), 
respectively,  and  are  expressed  as  the  sum  of  the  contributions  from  the  trunk  and  that  from 
the  branches. 


St(rt)=2%rJ+f2nrpl/!0rp 


(11) 


'min 


Vt(rt)=nr?l  +  /  v:rpl/f/lrp 


(12) 


Using  Eq.(10)  for  Nt,  St(rJ  and  Vt(r,)  become 
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(13) 


Sfr)=2vrft 


/ 

f  r,  Y| 

p 

t 

T  . 

\  mm// 

(14) 


where  the  (1-0)  term  in  St  has  been  included  to  account  for  pore  combination  (Simons, 
1979a). 


The  surface  area  associated  with  the  pore  tree  may  be  several  orders  of  magnitude 
greater  than  the  surface  area  of  the  trunk.  However,  the  volume  of  the  pore  tree  may,  at 
most,  be  one  order  of  magnitude  greater  than  that  of  the  trunk.  It  should  also  be  noted  that 
the  above  expressions  for  St  and  Vt  reduce  to  those  appropriate  to  a  single  cylindrical  pore  in 
the  limit  of  rt  ->rmin  (the  leaf  of  the  tree).  Furthermore,  the  integrals  of  St(rJ  and  V^rJ  over 
all  g(rt)  recover  the  total  internal  surface  area  and  pore  volume  of  the  porous  sample. 

Each  trunk  of  radius  rt  is  associated  with  a  specific  tree-like  structure  with  continuous 
branching  to  ever  decreasing  pore  radii.  The  radius  and  number  of  pores  is  a  unique  function 
of  the  distance  x  into  the  tree.  The  coordinate  x  is  skewed  in  that  it  follows  a  tortuous  path 
through  the  branches  of  the  tree.  Let  n(x)  represent  the  number  of  pores  of  radius  rp  at 
location  x  in  a  tree  of  trunk  radius  rt.  An  analysis  (Simons,  1982)  of  this  pore  tree  has 
demonstrated  that 


n(x)=r?/rfo) 


and  the  coordinate  x  is  related  to  rp  by 


drjdx=-rjlt 


(15) 


(16) 


The  continuous  branching  model  has  been  used  to  successfully  describe  char  oxidation 
(Lewis  and  Simons,  1979;  Simons,  1982  &  1983a),  coal  pyrolysis  (Simons,  1983b  &  1984) 
and  the  catalytic  cracking  of  benzene  by  porous  iron  oxides  (Simons  et  al.,  1986).  It  was  also 
used  to  successfully  describe  sulfur  sorption  (S02  and  H2S)  by  porous  calcine  (CaO)  in  the 
limit  of  zero  utilization  (Simons  and  Rawlins,  1980;  Simons  et  al.,  1984)  and  was  later 
extended  to  include  CaS04  and  CaS  deposits  (Simons  and  Garman,  1986;  Simons  et  al.,  1987; 
Simons,  1988;  Simons  et  al.,  1988).  The  subsequent  determination  of  the  controlling  physical 
parameters  led  to  a  new  concept  for  the  optimization  of  the  sulfur  sorption  process  (Simons, 
1991;  Simons  et  al.,  1992)  through  spray  drying  of  water  soluble  organic  calcium  solutions  to 
control  the  sorbent  pore  structure. 
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III.  INTERCONNECTIVITY  AND  PERMEABILITY 

The  first  step  in  determining  the  size  distribution  of  the  interconnected  pores  and  the 
distribution  of  the  permeability  is  to  determine  the  distribution  function  Gt(rp,r,)  drp  which 
represents  the  number  of  pores  of  radius  rp  (within  size  range  drp)  per  unit  cross  section  of  an 
arbitrary  plane  and  also  contained  within  a  tree  whose  trunk  radius  is  rt.  Consider  an  infinite 
homogeneous  isotropic  porous  media  and  isolate  a  spherical  volume  of  that  media  denoted  by 
the  radius  a.  Such  a  volume  is  illustrated  in  Fig.  2.  The  total  number  of  pores  of  radius  rp 
(within  size  range  drp)  intersecting  plane  AA  of  area  it  a2  has  previously  been  defined  by 
g(r^%a2dr  .  The  pores  in  plane  A  A  in  this  size  range  may  also  be  determined  by  integrating 

G(rt,rp)na2dr  over  all  trees  whose  trunk  intersects  the  exterior  surface  of  the  porous 
sample.  Hence  it  follows  that 


g(rp)na2drp=  f  [G(rt,rJ)na2drplg(rt)4na2drt 


where  only  those  trees  whose  trunk  radius  is  greater  than  rp  may  contain  a  pore  of  radius  rp . 

A  solution  to  Eq.  (17)  for  Gt(rp,r,)  will  not  necessarily  be  unique.  Physical  arguments 
will  help  determine  Gt(rp,r,)  and  help  ensure  that  it  is  the  particular  solution  we  seek.  Since  N, 
represents  the  number  of  pores  of  size  rp  in  the  tree  and  the  probability  of  a  pore  intersecting 
a  plane  is  proportional  to  its  length,  it  follows  that  G,(rp,r,)  should  be  proportional  to  the 
product  of  Nt  and  1/1, ,  i.e.,  proportional  to  r,2/rp3.  Eq.  (17)  is  identically  satisfied  by  a 


Figure  2.  Spherical  Volume  of  a  Porous  Media 


7 


function  which  differs  from  rt2/rp3  by  ln(rp). 


47ta2rjln(rmax/rp 


(18) 


Note  that  the  ln(rp/rmax)  introduces  an  integrable  singularity  at  rp=rmax.  No  physical 
significance  to  this  singularity  has  yet  been  identified. 

Further  confirmation  of  the  form  for  Gt(rp,r,)  may  be  made  through  the  definitions  of 
porosity  and  internal  surface  area.  If  the  porosity  and  internal  surface  area  in  plane  AA  due  to 
a  single  tree  of  radius  rt  are  defined  by 


and 


e,  -  / 


t,rp)2nrpdrp 


(19) 


(20) 


respectively,  then  the  total  porosity  and  internal  surface  area  must  be  represented  by  their 
respective  integrals  of  0t  and  Sj  over  all  trees.  Hence  0 


0  =  /  Qtg(rt)4na2dr( 

raia 


(21) 


and  sp 


sp  =  /  stg(rt)4na2drt 


'min 


(22) 


follow  respectively. 

It  is  stated  without  further  proof  that  Eqs.  (19)  to  (22)  are  compatible  in  both  their 
exact  and  their  approximate  forms.  The  distribution  function  Gt(rp,r,)  as  given  by  Eq.  (18)  is 
exact  and  will  be  used  to  determine  the  probability  of  trees  sharing  common  branches,  i.e., 
the  interconnectivity  of  the  pore  structure. 
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IV.  SUMMARY 


This  report  is  the  first  technical  report  on  the  project  entitled  "Development  of  Pore 
Structure  Models  for  Water  and  Contaminant  Transport  in  Partially  Frozen  Soils,"  submitted 
to  ARO  for  the  period  28  Sept  1994  to  13  October  1994.  The  pore  tree  model  is  being 
extended  to  describe  the  permeable  pore  structure  which  characterizes  the  subsurface  flow  of 
water  in  soil,  the  dispersion  of  contaminants  and  the  in-situ  remediation  of  contaminated  sites. 
Permeability  requires  a  statistical  determination  the  "branches"  that  are  common  to  several 
trees  to  allow  percolation  through  the  large  scale  (mobile)  structure  in  addition  to  diffusion 
through  the  smaller  scale  (immobile)  structure.  While  permeability  is  dominated  by  the  largest 
pores,  it  is  important  to  determine  the  level  of  convection  that  is  occurring  at  the  intermediate 
scales  in  order  to  accurately  describe  transport.  The  resulting  pore  structure  will  provide  an 
analytic  description  of  an  underground  network  upon  which  transport  and  coupled  chemical 
reactions  may  be  accurately  superimposed.  The  first  step  in  describing  the  permeability  is  to 
determine  the  pore  distribution  function  Gt(rp,r,)  drp  which  represents  the  number  of  pores  of 
radius  rp  (within  size  range  drp)  per  unit  cross  section  of  an  arbitrary  plane  and  also  contained 
within  a  specific  tree.  The  derivation  of  the  distribution  function  G,(rp,r,)  drp  is  described  in 
Section  III.  Subsequent  technical  reports  will  utilize  this  distribution  function  to  determine  the 
interconnectivity  of  the  pores  and  the  distribution  of  the  permeability  with  pore  size. 
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